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ABSTRACT 



Context. The X-ray transient MXB 0656-072 is a poorly studied member of high-mass X-ray binaries. Based on the transient nature 

of the X-ray emission, the detection of pulsations, and the early-type companion, it has been classified as a Be X-ray binary (Be/XRB). 

However, the flaring activity covering a large fraction of a giant outburst is somehow peculiar. 

Aims. Our goal is to investigate the multiwavelength variability of the high-mass X-ray binaiy MXB 0656-072. 

Methods. We carried out optical spectroscopy and analysed all RXTE archive data, performing a detailed X-ray-colour, spectral, and 

timing analysis of both normal (type-I) and giant (type-II) outbursts from MXB 0656-072. 

Results. This is the first detailed analysis of the optical counterpart in the classification region (4000-5000 A). From the strength 
and ratio of the elements and ions, we derive an 09.5Ve spectral type, in agreement with previous classification. This confirms its 
Be nature. The characterisation of the Be/XRB system relies on Balmer lines in emission in the optical spectra, long-term X-ray 
variability, and the orbital period vs. spin period and EW(Hq-) relation. The peculiar feature that distinguishes the type-ll outburst 
is flaring activity, which occurs during the whole outburst peak, before a smoother decay. We interpret it in terms of magneto- 
hydrodynamic instability. Colour and spectral analysis reveal a hardening of the spectrum as the flux increases. We explored the 
aperiodic X-ray variability of the system for the first time, finding a correlation of the central frequency and rms of the main timing 
component with luminosity, which extends up to a "saturation" flux of 1 x 10"** erg cm"^ s"' . A correlation between timing and spectral 
parameters was also found, pointing to an interconnection between the two physical regions responsible for both phenomenologies. 

Key words. X-rays: binaries — pulsars: individual: MXB 0656-072 



1. Introduction 

Be/X-ray binaries (Be/XRBs) constitute a sub-class of high- 
mass X-ray binaries (HMXBs) in which the companion is a 
Be star, i.e. a non-supergiant fast-rotating OB-star that during 
its life has shown at some point spectral lines in emission (see 
lReidl20l"Tl for a recent review). They are also characterised by 
infrared excess, which means that they are brighter in the IR 
than their non-emitting counterparts of the same spectral type. 
Both phenomena, emission lines and IR excess, are thought to 
arise from a common cause, namely the presence of an extended 
circumstellar envelope around the stellar equator, made up of 
ionised gas that is expelled from the star in a way that is not 
yet completely understo od. This complex scenario is referred to 
as the Be phenomenon (iPorter & Riviniusll2003t lEkstrom et al.l 
I2008h . 

When a Be star is part of an X-ray binary, the system is usu- 
ally transient, and the compact object is virtually always a pulsar, 
with typical spin periods ranging between 1-10^ s. Be/XRBs are 
characterised by high variability on a wide range of both time 
scales (from seconds to years) and wavelengths, although the 
fastest variability is observed in the X-ray band. For longer pe- 
riods, the variabiUty is apparent in both high-energy and low- 
energy wavelengths, and is attributed to major changes in the 
circumstellar disc structure. The complexity of the dynamics of 



the Be phenomenon and its relation with the accretion onto the 
compact object clearly require a multiwavelength approach in 
the study of these systems. 

Even if its phenomenology is entangled and multi-faceted, 
the long-term X-ray variability in Be/XRBs is traditionally de- 
scribed by a classification into two types of outbursts. Type-I (or 
normal) outbursts are periodic or quasi-periodic events, occur- 
ring in correspondence (or close) to the periastron passage of 
the neutron star They are generally short, with a typical dura- 
erg 



tion of 0.2-0.3 Poib, and show luminosities Lx - 10^^-10^^ 
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s" . Type-II (or giant) outbursts are unpredictable, long (one or 
more orbital periods), and bright events, with typical X-ray lu- 
minosities of Lx = lO^'^-lO^^ erg s"', coiTesponding to up to the 
Eddington luminosity for a neutron star The presence of quasi- 
periodic oscillations (QPOs) in some systems would support the 
suggestion of the formation of an accretion disc arou nd the neu- 
tron star during type-II outbursts (see for instance Motch et al.l 
llQQiHHavasaki & Okazakill2004 . 

The transient X-ray binary MXB 0656-072 was discov- 
ered by SAS-3 in September 1975, when a flux density of 80 
mCrab was reported (IClark et al.lll975l) . and subsequently ob- 
served twice in 1976 by A riel V at 50 (March 19) and 70 (March 
27) mCrab, respectively (iKaluzienskil 1 1 976h . These intensities 
would convert into an X-ray luminosity of ~2-3 x lO-'^ erg s"', 
assuming a distance of 3.9 kpc (McBiide et al..2006.) . Therefore, 
they would correspond to the typical X-ray luminosity range for 
type-I outbursts. 
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Table 1. Journal of RXTE observations. 



N. of 


Proposal 


MJD 


On-source 


pointings 


ID 


range 


time (ks) 


28 


80067 


52931.8-529754 


91.1 


33 


80430 


52966.9- 53033.3 


38.6 


44 


93032 


54419.5- 54748.6 


179.6 


123 


93423 


54449.1- 54776.6 


189.1 



Although the discovery of the source dates back to more than 
35 years ago, very little is known about the system. MXB 0656- 
072 was only catalogued as an HMXB in 2003 after extended 
re-brightening, when its optic al counterpart was identified and 
classified as an 09.7Ve star ("Pakull et al."l2003!), and a pulsed 
period of 160.7s detected (Morgan et al. 2003). The energy spec- 
trum of the source showed a cyclotron resonant ene rgy feature 
(CRSF) at a central energy of -33 keV (iHeindl etal.i200 3). This 
eve nt was classified as a t ype-II outburst. 

iMcBride et al.l (1200 6) monitored the source over the 2003 
major outburst and studied the average X-ray spectiTim during 
the peak of the outburst and the change in the spin period. They 
did a pulse phase-resolved analysis and found that the width 
of the CRSF varied with pulse phase, being wider during the 
pulse decline. However, the energy of the CRSF did not change 
with pulse phase. A cyclotron line at ~ 33 keV implies a mag- 
netic field strength of 3.7 x 10'^ G. They also present a timing 
analysis of the pulsations, including pulse profile dependence 
on energy and the changes i n the spin period throughout the 
outburst. Recentlv. I Yan et all (1201 2h have presented a correlated 
optical/X-ray analysis of the system during type-I outbursts, af- 
ter discovering a 101.2 days orbital period. 

In this paper we analyse all RXTE data available for MXB 
0656-072, which include the major 2003 outburst and a series 
of type-I outbursts observed in 2007-2008. We performed X-ray 
colour, spectral, and timing analysis of the giant outburst, and 
colour and spectral analysis of the normal outbursts. We focus on 
the aperiodic variability and the study of the broad-band noise. 
Moreover, we present here additional optical spectroscopy that 
allows us to refine and robustly justify the spectral classification 
on one side, besides providing a long-term correlated Ha/X-ray 
follow up of the system. 

2. RXTE observations and data reduction 

RXTE followed the source during the first giant outburst ob- 
served since its discovery, starting in October 2003 for approxi- 
mately three and a half months. (MJD 52931-53033). Renewed 
activity of the source was detected by the RXTE in November 
2007, lasting for one year (MJD 54419-54776), with luminosi- 
ties lower than in 2003. The 2007 X-ray variability consisted of 
a series of four (quasi-)periodic flares, which are reminiscent of 
type-I outbursts. Table 1 shows the observation log. The total net 
exposure amounted to 129.7 ks for the first event analysed here, 
and to 368.7 ks for the second one. 

We employed data fro m all the three instruments onboard 
RXTE dBradt etal.1 Il993h . the All-Sky Monitor (ASM), the 
Proportional Counter Array (PCA), and the High Energy X- 
ray Timing Experiment (HEXTE). The ASM incorporates three 
wide-angle shadow cameras equipped with proportional coun- 
ters with a total collecting area of 90 cm^. It works in the 2-10 
keV energy range, mapping the 80% of the sky every 90 min- 



utes. The PCA consists of five proportional counter units (PCUs) 
with a total collecting area of ~6250 cm^ and operates in the 2- 
60 keV range, with a nominal energy resolution of 18% at 6 keV. 
The HEXTE comprises two clusters of four Nal/Csl scintillation 
counters, with a total collecting area of 2 x 800 cm^, sensitive in 
the 15-250 keV band with a nominal energy resolution of 15% 
at 60 keV. Both the PCA and the HEXTE have a maximum time 
resolution of ~ I /is 

Data reduction was performed using HEASOFT version 6.9. 
An energy spectrum was obtained for each pointing, after filter- 
ing out unsuitable data according to the recommended critericQ 
employing Standard 2 mode data from the PCA (PCU2 only) 
and Standard (archive) mode from the HEXTE Cluster A, with 
a time resolution of 16s. The PCA and HEXTE spectra were 
extracted, background-subtracted, and dead-time corrected. For 
the PCA, the 3-30 keV energy range was retained, while the 
HEXTE provided a partially overlapping extension from 25 to 
100 keV. In the case of the 2003 outburst, for each observa- 
tion, the two resu lting spectra were simultaneously fitted with 
XSPEC V. 12.6.0 ( Amaudll9 96V For the 2007-2008 event, three 
average PCA spectra were extracted in three luminosity ranges 
(< 1 X lO-'', 1 - 2.5 X lO-'', and > 2.5 x 10"'^ erg cm^ s"' re- 
spectively) and fitted between 3-60 keV in order to constrain 
the CRSF. Then a spectiTim for each pointing was fitted between 
3-30 keV, fixing the CRSF parameters to the coiTesponding av- 
erage ones. HEXTE 2007-2008 spectra were too faint to be em- 
ployed. During the fitting, a systematic error of 0.6% was added 
to the PCA spectra. 

Power spectral density (PSD) was computed using PCA 
Event or Single_Bit data. We first extracted, for each observa- 
tion, a light curve in the energy range ~3.5-17 keV (channels 
8-39) with a time resolution of 2"^ s. The light curve was then 
divided into 128-s segments, and a fast Fourier transform was 
computed for each segment. The final PSD was computed as 
the average of all the power spectra obtained for each segment. 
These averaged power spectra were logarithmically rebinned in 
frequency and correcte d for dead time eff ects according to the 
prescriptions given in Nowa k et al.l (1 19991) . Power spectra were 
normalised such that the integral over the PSD is equal to the 
squared fractional rms amplitude, according to the so-called 
rmi^normalisation (iBelloni & Hasinge^ 1 19901: iMivamoto et al.l 
Il991h . 



3. Optical observations 

Optical spectroscopic observations of the companion star to 
MXB 0656-072 were performed using the Fred Lawrence 
Whipple Observatory at Mt. Hopkins (Arizona, USA) and the 
1.3-m telescope from the Skinakas observatory (Crete, Greece). 
Table|2]gives the log of the observations. 

The 1.3m telescope of the Skinakas Observatory (SKO) was 
equipped with a 2000 x 800 ISA SITe CCD and a 1302 1 
mm"' grating (on 30 September 2010) and 2400 1 mm"' (on 8 
November 201 1), giving a nominal dispersion of ~ 1 A/pixel and 
~ 0.5 A/pixel, respectively. We also observed MXB 0656-072 in 
queue mode with the 1.5-m telescope (FLWO) at Mt. Hopkins 
(Arizona) and the FAST-II spectrograph plus FAST3 CCD, a 
back-side illuminated 2688 x 512 UA STA520A chip with 15 
IJ.m pixels. Spectra of comparison lamps were taken before each 



' Among which, elevation from the Earth greater than 
10° and pointing offset lower than 0.02°; see PCA digest at 
http://heasarc.gsfc.nasa.gov/doc s/xte/pca-news.htmli 
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Table 2. Log of the optical observations. 
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JLf 

(2,400,000+) 


leiescope 


Orel ting 
(i/mm) 


Wavelength 
range (A) 


n vv \nct ) 
(A) 


tiw ynpj 
(A) 


14-11-2009* 


55150.41 


FLWO 


600 


4760-6760 


-18.4+1.0 


-2.98 + 0.11 


12-01-2010 


55209.37 


FLWO 


600 


4730-6730 


-20.9+1.5 


-3.43 + 0.08 


16-01-2010 


55213.21 


FLWO 


600 


4740-6740 


-21.2 + 0.8 


-3.60 + 0.09 


30-09-2010 


55470.60 


SKO 


1301 


5300-7300 


-10.8 + 0.4 


- 


30-10-2010 


55501.01 


FLWO 


1200 


6200-7200 


-11.9 + 0.6 


- 


z9-i i-zUiU 


JJOiU.oi 


rLWU 


1 OAn 


ozUO- //OU 


1 1 1 A /C 

-t 1.0 + 0.0 




03-10-2011 


55838.98 


FLWO 


1200 


6200-7200 


-12.2 + 0.7 




01-11-2011 


55867.88 


FLWO 


1200 


6200-7200 


-13.0 + 0.7 




08-11-2011 


55873.60 


SKO 


2400 


3940-5040 




-2.23 + 0.07 


23-11-2011 


55889.93 


FLWO 


1200 


6200-7200 


-14.1+0.8 




31-11-2011* 


55927.71 


FLWO 


1200 


6200-7200 


-15.4 + 0.8 




19-01-2012 


55946.79 


FLWO 


1200 


6200-7200 


-15.3+0.8 




22-01-2012* 


55927.71 


FLWO 


1200 


6200-7200 


-15.2 + 0.8 





*: Average of two measurements. 



exposure to account for small variations in the wavelength cali- 
bration during the night. To ensure a homogeneous processing of 
the spectra, all of them were normalised with respect to the local 
continuum, which was rectified to unity by employing a spline 
fit. 

We measured calibrated photometry of the optical counter- 
part with dedicated observations for the first time, performed 
from the 1.3-m telescope of the Skinakas Observatory on 2 
November 2010 (JD 2,455,503.6). MXB 0656-072 was observed 
through the Johnson B, V, R, and I filters. The telescope was 
equipped with a 2048 x 2048 ANDOR CCD with a 13.5 yum 
pixel size. Standard stars from the Landolt list (Landolt 2009) 
were used for the transformation equations. Reduction of the 
data was carried out in the standard way using the IRAF tools 
for aperture photometry. The resulting magnitudes are: B - 
13.25 + 0.02 mag, V = 12.25 + 0.02 mag, R = 11.63 + 0.02 
mag, and / = 10.97 + 0.02 mag. 

4. Spectral class 

The only report of the spect ral type of th e mass ive companion in 
MXB 0656-072 is given bv iPakull etal](l2003l) . who suggest an 
09.7V spectral class. Blue- and red-end op tical spec tra covering 
the period 2005-2009 are also presented in lYan et al.l(l2012l) . 

Figure [T] shows the optical spectrum of MXB 0656-072 in 
the region 4000-4800 A from the Skinakas observatory. This 
spectrum is the average of three spectra obtained with a total 
exposure time of 7200s each. The distinct presence of He II 
lines indicates an O-type star, while the presence of He I lines 
implies that the spectral type must be later than 08. Some He 
I lines, such as A4113 and A4921, appear (partially) in emis- 
sion and cannot be separated out from the continuum. The ra- 
tio C III A 4650 t o He II A46S6 is close to 1, w hich agrees with 
an 09-09.5 type (IWalborn & FitzpatrickllT990l) . The ratio He II 
/14200/He I /i4144 allows us to distinguish between these two 
close sub-classes (Walborn 1971). This ratio is approximately 1 
in 09 stars and lower than 1 in 09.5 stars. In MXB 0656-072 this 
ratio appears to be slightly lower than 1, favouring the later type 
classification. Also, since the strength of the He I /14144 might 
be diminished by the emission from the circumstellar disc, the 
09.5 class appears to be more likely. On the other hand, given 
the relatively low S/N in this part of the spectrum and the uncer- 
tainty introduced in the definition of the continuum during the 



normalisation, an 09 spectral type cannot be completely ruled 
out. 

As for the luminosity class, the strong He II /14686 absorp- 
tion accompanied by weak N III /14634-4640-4642 clearly in- 
dicates a main-sequence star, and so does the strength of Si IV 
/14089 in comparison with that of He I A4144. We conclude that 
the optical counterpart to MXB 0656-072 is an 09.5V stafl 

In addition to He I lines, the Balmer series of hydrogen lines 
are strongly affected by emission. Even He appears to be filled- 
in with emission. Columns 6 and 7 in Table |2] give the equiva- 
lent width of the Ha and li/3 lines. A long-term decrease in the 
strength of these two lines is observed. 

5. Results 

We analysed all the observations of MXB 0656-072 in the RXTE 
archive. We separated the observations into two intervals corre- 
sponding to two significant events. The first interval started on 
20 October 2003 (MJD 52932) and covered a total of 101.5 days. 
This interval incl udes a giant (type II) outburst. Assuming a dis- 
tance of 3.9 kpc dMcBride et al.ll2006l) . the maximum 3—30 keV 
luminosity during the outburst was Lx - 3.7 x lO-'^ erg s"', 
registered at MJD~52966.9. The second interval covered the pe- 
riod 27 November 2007 to 11 November 2008 (MJD 54419.5- 
54776.6) and includes a series of minor (type I) outbursts. The 
peak X-ray luminosity of these outbursts was Lx - 1.37 x 10^^ 
erg s"'. 

5. 1 . Type II outburst 
5.1.1. Colour analysis 

The PCA light curve and colour behaviour during the 2003 out- 
burst is presented in Fig. |2] The outburst showed strong flaring 
behaviour during the peak phase, followed by smoother decay. 
Each point corresponds to an RXTE pointing and is directly ob- 
tained from PCU2 count rate. Different symbols mark the differ- 

- The interpolated class 09.7 suggested by Pakull et alj ( l2003h is 
mainly used for supergiants. The prim ary definin g criterion for the 
09.7 type is He II ^4541^1 III ^4552 jWalbom & Fitzpatrick 1996|). 
According to W albornl(ll97lh . it does not appear that this subdivision is 
useful at the lower luminosities (class II and below) because the Si III 
lines are too weak. 
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Fig. 1. The optical spectrum of the optical counterpart to MXB 
0656-072 in the 4000-4800 A region. The identified lines corre- 
spond to He II ^4200, ^4541, A46E6, He I A4026, A4144, /14387, 
i4471, ^4713, C III ^4070, /14650, N III /14640, Si IV ^4089, 
/141 16, and the hydrogen lines of the Balmer series between H/j 
and Hf. A Gaussian smoothing filter (cr = 1) was applied to 
reduce the noise. 
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Fig. 2. PCA light curve and colour behaviour during the 2003 
outburst. In the inset, 2s resolution light curve for one pointing 
during the flaring phase. 



ent outburst phases, the flare-like phase (open squares), and the 
smoother decay (filled circles). Error bars are the same size as 
the points whenever they do not appear in the plots. The X-ray 
colours were defined as follows, soft colour (SC): 7-10 keV / 
4-7 keV; hai'd colour (HC): 15-30 keV / 10-15 keV. The two 
colours follow identical patterns during the outburst, both cor- 
relating with flux. In the inset, a zoomed-in view is shown, pre- 
senting the light curve for one PCA pointing, with a 2s time 



resolution. From the inset it is clear that the flare-like activity 
displayed during the outburst is large-scale behaviour that in fact 
corresponds to variability on various time scales, if investigated 
at higher time resolution: besides the ~160s pulse period, slower 
and faster changes in intensity are clearly detected in the 2s res- 
olution light curve. 

The amplitude of change in count rate in both colours is 
more than twice in the decay (~0. 15) compared to the flare phase 
(~0.05, see Fig.|2|. Also, the colour values are higher larger dur- 
ing the flares, indicating a harder spectrum. 

5.1.2. Spectral analysis 

We fitted energy spectra with a continuum constituted by a 
photo-absorbed power law with a high-energy exponential cut- 
off, modified by a Gaussian line at ~6.5 keV with a fixed 0.5 keV 
width to account for Fe Ko' fluorescence. A CRSF at an average 
central energy of 37 keV was detected in high-flux observations 
only, above Lx = 1-6 x lO-'^ erg s"', i.e. above 0.4 x Lxmax- 

In addition, MXB 0656-072 shows significant residuals 
at ~11 keV. They were fitted out by means of a Gaussian 
absorption-like profile, which allowed acceptable fits, passing 
from x'- -107 for 70 DOF to x'' ~85 for 67 DOF, for a 
typical high-flux observation. This component is found at al- 
most constant energy across the spectra, with a weighted mean 
value of 1 1.68±0.05 keV, and was consistently reported also by 
[McBride et al. (2006) for the giant outburst and by Yan et al] 
(l2012h for the normal ones. We found that this feature is only 
necessary in spectra where the CRSF is present as well. Its ori- 
gin is uncertain. 

Figure [3] shows a typical spectrum for a high-flux observa- 
tion, with the corresponding residuals in the case of the best fit 
(a), the fit excluding the absorption at ~ 1 1 keV (b), and flie one 
excluding the CRSF (c) from the model. 

The best-fit main spectral parameters are shown in Fig. 
as a function of the calculated 3-30 keV flux. Different marks 
and colours identify the two different phases of the outburst, the 
flaring and the smooth one. Points separated in time and corre- 
sponding to different phases behave in a similar way, which only 
depends on flux. The power-law photon index decreases with X- 
ray flux, confirming the result from the colour analysis that as 
the X-ray flux increases the spectrum becomes harder (see also 
Fig. |2]i. The strength of the iron fluorescence line nicely corre- 
lates with flux, revealing an increase in the reprocessed material 
as the flux increases, as expected. The central energy of the Fe 
line remained fairly constant, with a mean value of 6.49+0.06 
keV, and a mean equivalent width (EW) of 0.34±0.07 keV. 

The CRSF was almost constant during the outburst, with 
the following weighted mean values: = 36.8 + 0.4 keV, 
cr = 9.1 + .4. These value s are n ot compatible with the best- 
fit values bv lMcBride etalJ (l2006l) . who found the CRSF at a 
cent ral energy of 32 .8^°^ keV. The discrepancy may be because 



iMcBride et al. l (|2006) obtained one spectral fit from a spectrum 
retrieved by summing up all the spectra from the flare phase of 
the outburst (from MJD 52932 to MJD 52964). In fact, when 
choosing four individual observations of that interval, the cy- 
clotron line in the spectra of those four obs ervations s howed line 
energies above 33 keV (see Fig. 3 in McBride et al.il2006) . Our 
best-fit val ue for the central en ergy agrees with the one firstly 
reported bv lHeindl etan(l2003h . of = 36 + 1 keV. 

All the observed correlations, except for the iron line 
strength, are significant up to a "saturation" flux of ~10"^ erg 
cm"^ s"', which corresponds to Lx = 1.8 x lO-'^ erg s"^ This 
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Fig. 3. Typical spectrum for an observation during the 
flaring phase of the giant outburst (obsid: 80067-11-03-05, 
Lx=2.8xl0" erg s '): in the upper panel, PC A and HEXTE data 
points and corresponding best fit are reported; below: residuals 
for the best fit (a), the fit excluding the absorption feature at ~ 1 1 
keV (b), and the one excluding the CRSF (c), respectively. 
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Fig. 4. Evolution of the main spectral parameters during the 2003 
outburst. From the upper panel: photon index, cutoff energy, hy- 
drogen column density, and Fe line intensity. 



luminosity roughly coincides with the large-amplitude flaring 
phase. 

Because of the flaring activity, it is difficult to identify the 
peak of the outburst. The average 3-30 keV X-ray luminosity in 
the time interval MJD 52940-52970 is Lx = 2.5 x lO" erg s"', 
although a maximum of Lx - 3.7 x 10^^ erg s"' was obtained 
on MJD 52966.9. The lowest luminosity corresponds to the last 
point of the decay with Lx = 1x10^* erg s ' . 



5.1.3. Timing analysis 

In this work, we focus on the aperiodic variability o f the system . 
For a study of the X-ray pulsations, see McBride et al.l (l2006h . 
The neutron star spin frequency has a fundamental peak at ~6 
mHz, below our PSD frequency range (0.008-32 Hz). Thus, 
peaks derived from the neutron star's pulsations do not appre- 
ciably distort the continuum in the power spectra. 

We fitted each PSD with the sum of Lorentzian functions 
with the objective of providing a unified phenomenological de- 
scription of the timing behaviour of the system during the out- 
burst. We denote each component as Li, and its characteristic 
freque ncy Vmax as Vi. According to the definition in lBelloni et al.l 
(I2002h . this is the frequency where the component contributes 
most of its variance per logarithmic frequency interval, Vn,ax - 



1 - 
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+ {FWHM/2)^, where vo is the centroid frequency and 

FWHM is the full width at half maximum of the Lorentzian 
function. In this work, we always refer to characteristic frequen- 
cies Vmax- 

The low and middle-frequency noise (Lq and Li) is ac- 
counted for by zero-centred Lorentzians with a characteristic fre- 
quency that is generally lower for low-flux pointings, and higher 
for high-flux observations. The corresponding fractional rms 



* 











60 



m 40- 



••• 



10" 



3-30 keV flux (erg cm"^ s"^) 



10" 



Fig. 5. Evolution of the timing parameters, characteristic fre- 
quency and fractional rms, of the best constrained noise com- 
ponent, Li, during the 2003 outburst. 
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Fig. 6. Relation between the maximum frequencies of the Li and 
La components. 
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Fig. 7. Relation between the photon index and the characteristic 
frequency of the Li broad-noise component. 
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varies during the outburst, in anti-correlation with flux. These 
components are the only ones required during all the smooth de- 
cay phase (after MJD 52975) and part of the flaring phase of the 
outburst, up to /i a; 1 X 10"'^ erg cm"^ s"'. Beyond that luminos- 
ity, an additional component is necessary, L2, whose characteris- 
tic frequency varies in the range ~l-5 Hz, whereas the fractional 
amplitude of variability varies between 10% and 20%, without a 
clear dependence on flux. Unlike Lq and L\ , this component has, 
in general, a non-zero centroid frequency, and its average value 
for the g-factor, defined as g = vqIFWHM, is ~0.6, denot- 
ing a narrower feature compared to Lq and L\. Figure |5]presents 
the evolution of the characteristic frequency and rms over the 
outburst for the Li component, the best constrained one. Up to 

«i 1 X 10"^ erg cm"^ s"', as the flux increases the charac- 
teristic frequency increases, whilst after that luminosity it re- 
mains approximately constant. The rms variability shows an op- 
posite trend, decreasing as the flux increases, and then saturating. 
Figure|6]shows the relation between the L\ and L2 characteristic 
frequencies; although within some scattering, the two frequen- 
cies follow a correlated trend in all their range of variation. 

We studied the correlated spectral/timing behaviour, and 
found that during the very last part of the flaring phase, and the 
whole decay, the photon index F and the central frequency of 
the main timing component, vi, vary in an anti-correlated way 
(Fig. |7|. During the flares the L\ frequency shows only slight 
variation due to the appearance of L2, so that no spectral/timing 
relation could be expected. 



Fig. 8. ASM light curve and Ha equivalent width evolution. A 5- 
day running average was applied to the 1-day ASM light curve 
to reduce the noise . Open circles correspond to observations 
bv lYanetal.1 (l2012h and filled circles to our observations (see 
Table|2]l. 
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Fig. 9. f 01b - EW{Ha) diagram. The star symbol marks the po- 
sition of MXB 0656-072. 



5.2. Type I outbursts 

In addition to the major (type II) outburst reported in previous 
sections, MXB 0656-072 underwent a series of four fainter out- 
bursts between November 2007 and November 2008. The PC A 
began to monitor these outbursts at the end of the first one. All 
outbursts exhibited similar peak luminosities (Lpeak = 1.2 x lO-'^ 
erg s"'). This luminosity is about three times lower than that of 
the November 2003 type-II outburst. 

We show in Fig. |8] the simultaneous Ha measurements and 
ASM light curve. The outbursts are separated by ~100 days and 
in between the major peaks, other minor peaks are observed. 
The onset of these outbursts coincided with an optical maximum 
brightness of the donor: at around MJD 54500 (February 2008) 
the equivalent width of the Ha line seems to have reached a max- 
imum value of -25 A dYan et alJl20T2l) . 

In our X-ray spectral analysis, in order to best constrain the 
CRSF, hardly detectable at very low flux, we extracted three 
average spectra in three ranges of luminosities and fitted them 
between 3-60 keV. The model employed was the same as for 



type-II outburst, and all the components were needed to obtain 
acceptable fits at in all the three spectra. Spectral parameters are 
shown in Table |3] Similarly to our analysis of the giant outburst, 
we also generated one energy spectrum for each observing in- 
terval and studied the evolution of the spectral parameters . We 
fixed the CRSF parameters in each single spectrum to the ones 
obtained from the corresponding flux average spectrum, and fit- 
ted each one between 3 keV and 30 keV. We found that the spec- 
tral parameters follow similar trends to those seen during the 
type II outburst (Fig. [Hi. The photon index decreases as the flux 
increases, as in the type II outburst, although with a steeper de- 
pendence with flux. The iron line energy does not significantly 
vary, while its intensity, expressed as normalisation, tightly cor- 
relates with luminosity. The iron line width was fixed at 0.5 keV. 
The hydrogen column density generally anti-correlates with flux, 
while the cutoff energy does not show any smooth relation with 
luminosity, although it displays lower values at higher flux and 
vice versa. 
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Table 3. Type-I outburst spectral analysis for average spectra at 
different flux ranges. 



Spectral parameter low flux" med. flux* high flux'^ 



r 


1.08± 


:0.08 


83+0 04 


0.49±0.03 


cutoff en. (keV) 


16.9 


,+2.5 
-2.3 




10.5+0.02 


pow. norm. (ph/keV/cm^/s) 


o.ii: 


t0.04 
-0.02 


0.12±0.01 


115+'"'o-* 

-0.007 


nH (lO-^cm--) 


4.1±0.5 


2.8±0.3 


1.9+0.3 


Epe (keV) 


6.5i 


:0.1 




6.5+0.4 


EWpe 


0.06± 


:0.02 




0.17!!!:?? 


ElOkeV-gaus 


10.3 


+0.4 
-0.5 


10.6!°:] 


10.6!°:? 


O" lOkeV-gaus 


6 5+° '' 


4 7+0.3 
^•^-0.5 


9 Q+0.24 
^^-0.23 


TlOkcV-gaus 


15.8 


+6.2 
-3.9 


9 Q+0.4 
^•■'-0.5 


1 05+0'* 


F 


31.9 


+ 1.3 
-2.7 


35.6!?-^ 


35.0!!:^ 






8.5!!:^ 






58.0: 


+ 13.2 
-8.9 


14.0!^:^ 


13.0!^:^ 



a: < 10 ' erg cm- s ' 

b: (1-2.5) xlO-' erg cm- s"' 

c: > 2.5 X 10"' erg cm^ s"' 



6. Discussion 

We have performed a detailed X-ray and optical analysis of the 
poorly studied hard X-ray transient MXB 0656-072. All the 
available observational data indicate that MXB 0656-072 is a 
member of the class of massive X-ray binaries known as Be/X- 
ray binaries. X-rays are produced in the vicinity of the compact 
object, while the optical variability comes from the young and 
massive companion. The detection of X-ray pulsations, the tran- 
sient nature of the X-ray emission, and the characteristics of 
the X-ray spectrum (power-law continuum modified by an ex- 
ponential cutoff and the presence of fluorescence iron line and 
cyclotron feature) are typical of neutron star binaries. The obser- 
vation of Balmer lines in emission favours the Be/XRB classifi- 
cation. Furthermore, the long-term X-ray variability, consisting 
of giant (type II) and minor recurrent outbursts (type I) are typi- 
cal of Be/XRB. 

6.1. Optical observations 

The optical counterpart to MXB 0656-072 was classified as 
a 09.5 Ve star, refining previous classification by iPakuU et"al] 
(l2003h . Its spectrum is strongly affected by emission, with the 
first three lines of the Balmer series (Hot, Hj6, and Hy) showing 
an emission profile, while the next two (H5 and He) are par- 
tially filled in with emission. This extra emission is thought to 
arise from the equatorial disc around the Be star. The picture 
described by optical spectroscopy is fully consistent with the 
so-called "Be-phenomenon" and confirms that the system is a 

Be/XRB. 

From the analysis of type-I outbursts, lYan et al] (1201 2l) find 
an orbital period of 101. 2d. Once the period is known, we can 
use two important relationshi ps inv olving the orbital period of 
the system thePsmn - f orb (iCorbelill98Q) and Porb - EW{Ha) 
(iReig et al.ll997t lR ei j201lh diagrams, to s upport the orbita l pe- 
riod found. In the first diagram (see Fig. 6 in lYan et al.l2012l) . the 
source is clearly located in the region occupied by Be/XRBs; in 
the second one (Fig.|9]l, a ~100d orbital period fits nicely in the 
expected EW(Hq') vs. orbital period relation. The Po^-EW{Ha) 
correlation is a consequence of tidal truncation of the Be star's 



circumstellar disc by the neutron star. Assuming that the equiva- 
lent width of the Ha line, EW{Ha), provides a good measure 
of the size of the circumstellar disc (Ouirrenbach et al. 199^ 
iTvcner eF al. 2005), the Poj\,-EW{Ha) correlation indicates that 
systems with long orbital periods have larger discs, while narrow 
orbit systems contain smaller discs. In short orbital period sys- 
tems, the neutron star prevents the formation of extended discs. 
As can be seen in Fig.|9l an orbital period of ~ 100 days agrees 
very well with the maximum EWiHa) of 25 A. 

6.2. Type-ll outburst 

The profile of the type-II outburst in MXB 0656-072 is unusual 
among Be/XRB. Although flaring behaviour h as been seen in 
other systems (e.^., EXO 2() 30-h375, A0535-H26. lKlochkov et al] 
1201 UlCaballero et al.ll2008l) . the profile of the outbursts tends to 
be more symmetric and to have a smoother peak. In MXB 0656- 
072, the flaring activity covers a substantial part of the outburst 
(Fig. |2]l. It was found that magneto-hydrodynamic instabilities 
at the inner edge of the accretion disc may produce oscillations 
in the accreting flow, possibly leading to the observed behaviour 
(Apparao 1991; Postnov et al. 2008; D'Angelo & Spruit 2010|. 
Especially, D'Angelo & Spr uitI (120101) show that, if the accretion 
disc is truncated by the neutron star's magnetic field outside but 
close to the corotation radius (at which the Keplerian frequency 
in the disc equals the star's rotational frequency), then the ac- 
cretion becomes time-dependent and takes the form of repeated 
bursts. 

We found correlated behaviour of both the SC and the HC 
with flux, which corresponds to a general hardening of the spec- 
tra as the flux increases. The dependence of the SC on luminosity 
does, in appearance, contrast with th e observed negative corre- 
lation in other sources. iReigl (l2008l) analysed the spectral and 
timing properties of four Be/XRBs during giant outburst, iden- 
tifying two source states, a low-luminosity horizontal branch 
(HB) and a high-luminosity diagonal branch (DB). During the 
outburst, the sources spend most of the time in the DB, during 
which the SC anti-correlates with flux. The transition to the HB 
inverts this trend, and the SC starts correlates with it. The corre- 
lation between the SC and luminosity observed in MXB 0656- 
072 thus reveals that the source never undergoes the transition 
to the DB, but remains in the HB during all the duration of the 
outburst. 

The hardening of the spectrum at high flux is confirmed by 
spectral analysis: from simultaneous spectral fitting of PCA and 
HEXTE spectra, we found a clear anti-correlated behaviour of 
the spectral index with flux (Fig. lU up to flux ~ 1.2 x 10"^ 
erg cm^ s"', or Lx = 2.3 x 10^^ erg s"', when the trend flat- 
tens. The spectral hardening with increasing flux can be inter- 
preted with simple Comptonisation models: the X-ray emission 
from HMXBs is thought to arise in an accretion column, close 
to the neutron star's surface, as a result of Comptonisation of 
soft photons injected in the accretion flow from the NS ther- 
mal jnoimdj_bj;_Jiigl>energy electrons of the accreting mat- 
ter ( Becker & WoTllIOOTh . Increasing X-ray luminosity would 
mainly correspond to increasing mass accretion rate, which can 
be translated into averagely increasing number of up-scattering 
collisions of photons with electrons, resulting in a harder spec- 
trum. 

After lA 1118-615 (iNespoU & Rei^ll201 ih . MXB 0656-072 
is the second X-ray pulsar to show correlated spectral/aperiodic 
behaviour (Fig|7]i. This translates into a necessary connection 
between the two physical regions responsible for producing the 
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two phenomenologies, the accretion column on one side, where 
energy spectra arise, and the accretion disc on the other side 
where, according to the "perturbation propagatio n" model, the 
aperiodic variability is located (iRevnivtsev et al.ll2()09i and ref- 
erences therein). In this model, in fact, the X-ray variability is 
caused by perturbations in the inner disc flow, at different radii. 
The two regions are physically separated, being the first close to 
the NS polar cap regions, and the second confined outside the 
magnetosphere, but our results on two systems show that they 
are somehow coupled. 

In general, the X-ray variability of MXB 0656-072 during 
the type-II outburst resembles that of 1 A 1 118-615 at many lev- 
els. The colours behave exactly the same in the two sources (both 
correlate with flux). The spectral parameters have the same trend 
in the two sources, and in both cases the relation with flux ceases 
at some saturation luminosity. In MXB 0656-072 this saturation 
is reached at 0.5xLxmax, while in lA 1118-615 at 0.7xLxmax- 
Finally, the correlations between the two broadband timing com- 
ponents and between spectral and timing parameters are also 
observed in both sources, making them very similar during a 
type-II outburst, although no flaring activity is observed in lA 
1118-615. 

A peculiar feature in the X-ray energy spectra of MXB 0656- 
072 is an absorption-line-like profile at an average constant en- 
ergy between 11-12 keV. This component was investigated by 
[Coburn (2001), who detected it in the range 8-12 keV in the 
spectra of many X-ray pulsars. The feature was consistently ob- 
served at the same energy, irrespective the CRSF energy, and 
moreover, it was evident in some systems that do not display a 
cyclotron line. This made Coburiji (120011) conclude that the com- 
ponent should not be a magnetic effect. The feature seems to 
be intrinsic to X-ray pulsars spectra, since it was observed with 
different instrume nts (besides RXTE , Ginga and BeppoSAX, 
[Mihara 1995; Sant angelo et al.ll998l) . In the case of MXB 0656- 
072, this component is only found in spectra where the CRSF is 
present as well, although no other relation could be established 
between the two features. 

7. Conclusions 

We presented a detailed X-ray and optical study of MXB 0656- 
072 covering both types of X-ray variability observed in a 
Be/XRB, namely type-I and type-II outbursts. The major out- 
burst is characterised by flare-like behaviour during its pe ak, fol- 
lowed by smoother decay. We interpreted the flaring activity as 
possibly due to magneto-hydrodynamic instabilities at the inner 
edge of the accretion disc. The colour and spectral analyses re- 
veal a hardening of the spectra as the luminosity increases, which 
can be understood in the framework of the models for spectra 
production in X-ray pulsars. The analysis of aperiodic variability 
shows correlated behaviour of the timing parameters with flux, 
which translates into a correlation between spectral/timing fea- 
tures and can be interpreted as an interconnection between the 
two physical regions responsible for the two phenomenologies. 
All the X-ray behaviour during the type-II outburst resembles 
that of lA 1118-615, although no such flaring activity was ob- 
served in that system. The spin period vs. EW(Ha) relation con- 
firmed the orbital period proposed for the source. The full multi- 
wavelength analysis corroborates the Be/XRB nature of the sys- 
tem. 

Further observations during major outbursts are needed in 
order to explore the nature of the flaring activity and the tim- 
ing/spectral correlation detected in this work. Deeper compre- 
hension of the interaction between the magnetosphere and the 



accretion disc is also necessary to explain the correlated be- 
haviour of the spectral and aperiodic features. 
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